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Abstract 

We study higher dimensional models with magnetic fluxes, which can be derived from 
superstring theory. We study mass spectrum and wavefunctions of massless and mas- 
sive modes for spinor, scalar and vector fields. We compute the 3-point couplings and 
higher order couplings among massless modes and massive modes in 4D low-energy ef- 
fective field theory. These couplings have non-trivial behaviors, because wavefunctions 
of massless and massive modes are non-trivial. 



1 Introduction 



Field theory in higher dimensions plays a role in particle physics and cosmology. In partic- 
ular, extra dimensional field theory derived from superstring theory is important. Their 
four-dimensional (4D) low-energy effective field theories are determined by geometrical 
aspects of compact extra dimensions. One of the simplest compact spaces is a torus. 
However, the simple toroidal compactification does not lead to a chiral theory as a 4D 
low-energy effective field theory. Hence, it is a key issue to realize a 4D chiral theory when 
we start with higher dimensional field theory. 

Complicated geometrical backgrounds such as Calabi-Yau manifolds would lead to a 
4D chiral theory, although it may be difficult to compute explicitly 4D low-energy effec- 
tive field theories from such geometrical backgrounds. On the other hand, the toroidal 
compactification can also lead to a 4D chiral theory when we introduce non-vanishing 
magnetic fiuxes in extra dimensions. The number of zero-modes are determined by the 
size of magnetic fiux and each zero-mode has a quasi-localized profile. Thus, the toroidal 
compactification with magnetic fiuxes is quite attractive background for higher dimen- 
sional field theory [U El El IH IS] (see also [SIE])- Its stringy setup corresponds to magnetized 
D-brane models wrapping cycles on the torus [SI El [IDl E] • Furthermore, magnetized D- 
brane models are the T-dual of intersecting D-brane models, and many interesting models 
have been constructed in both types of models [121 113] . 

The Yukawa couplings among massless modes were computed by integrating the over- 
lap of wavefunctions in the extra dimensional space [1] . If zero-modes are quasi-legalized 
far away from each other, their couplings are suppressed. On the other hand, if they 
are localized near each other, their couplings are not suppressed, but would be of C(l). 
Thus, these localization behaviors are important from the phenomenological viewpoint, 
for example, to derive the realistic values of quark and lepton masses and their mixing 
angles. Furthermore, higher order couplings among massless modes were also computed 
|14] . Interestingly, they are written by products of 3-point couplings. These low-energy 
effective field theories can also lead to Abelian and non-Abelian discrete fiavor symme- 
tries, e.g. and A(27) fiavor symmetries [T5l [TBjFI These non-Abelian discrete fiavor 
symmetries are important to derive the realistic quark and lepton mass matrices (see e.g. 
[18] and references therein). 

In addition to massless modes, massive modes also have important effects in 4D low- 
energy effective field theory. For example, they may induce the fast proton decay and 
fiavor changing neutral currents (FCNCs) (see e.g. PJIJ). Our purpose in this paper 
is to study massive modes in the extra dimensional models with magnetic fiuxes. We 
study their mass spectrum and wavefunctions explicitly. Then, we study compute 3-point 
couplings and higher order couplings including these massive modes. These couplings 
have non-trivial behaviors, because wavefunctions of massless and massive modes are 
non-trivial. 

This paper is organized as follows. In section 2, we briefiy review on the fermion zero- 
modes on with the magnetic fiux. Then, we study mass spectrum and wavefunctions 

^ Similar non-Abelian discrete flavor symmetries are derived in heterotic orbifold models [17j . 
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of higher modes exphcitly. These analysis is extended to those for zero-modes and higher 
modes of scalar and vector fields. Its extension to is straightforward. In section 3, 
we compute couplings among these modes. In section 13. we give a brief review on 
computations of the 3-point couplings and higher order couplings among zero-modes. 
Then, we extend them to the computations of the 3-point and higher order couplings 
including higher modes in 13.21 In section 13. 3[ we also consider the couplings including 
massive modes due to only the Wilson line effect, but not magnetic fluxes. In section 
4, we give comments on some phenomenological implications of our results. Section 5 is 
devoted to conclusion and discussion. In Appendix A, we show some useful properties 
of the Hermite function. In Appendix B, we briefly review on the vector field in extra 
dimensions. In Appendix C, we show useful properties of the products of zero-mode 
wavef unctions. 

2 Mass spectrum and wavefunctions of massive modes 

We consider the (4 + c?)-dimensions, and denote four-dimensional and (i- dimensional co- 
ordinates by x'^ and with /i = 0, ■ ■ ■ ,3 and m = 1, ■ ■ ■ ,d, respectively. We study 
the spinor field X{x^,y'^) and the vector field AM{x^,y"^) with M = 0, ■ ■ ■ , (3 + (i). We 
decompose these fields as follows, 

A(x^n = J^^nix^Myn^ (2.1) 

n 

AM{x^,yn = 5^<^n,M(x^)0n,M(l/™). (2.2) 
n 

Here we choose the internal wavefunctions ipn{y"^) as eigenfunctions of the internal Dirac 
operator as 

iT'^D^ij^ = m,,ij„, (2.3) 

where F*" denote the gamma matrices in the internal space. The eigenvalues of m„ become 
masses of the modes Xn{x'^) in 4D effective field theory. Similarly, 4>n,M{y"^) correspond 
to eigenfunctions of the internal Laplace operators, as will be shown explicitly later. The 
scalar field in the (4 + (i)- dimensions is also decomposed in a similar way. 

2.1 with magnetic flux 

First, let us consider the 2D torus, T^. Here, we follow the notation of Ref. [1]. Instead of 
the real coordinates y^ and y"^, we use the complex coordinate, z = y^ + ry'^ with r G C. 
The metric is given by 

ds^ = 2{27iRfdzdz. (2.4) 

We identify the complex coordinate as z ^ z + 1 and 2; ~ z + r on T^. The area is written 
by ^ = An'^RHmT. 
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We introduce the U(l) magnetic flux on as 



m. (2.5) 



imr 

This magnetic flux is derived e.g. from tlie following vector potential, 

Az = — mz, Az = — mz. (2.6) 

2Imr ' 2Imr ^ ^ 

Their boundary conditions can be written as 

A{z + 1) = Ai{z) + diXi, A{z + t) = A{z) + d.a2, (2.7) 

where 

71 TT 

Xi = m Imz, X2 = ITT' Imf^;. (2-8) 

Imr Imr 

Furthermore, we can introduce non-vanishing Wilson lines by using 

TT TT _ 

Xi = Im(m2; + a), X2 = lmf{mz + a), (2-9) 

Imr Imr 

where a is complex and corresponds to the degree of freedom of the Wilson line. It is 
convenient to use the following notation, 

a = mC, (2.10) 

for m 7^ 0. 

2.1.1 Fermion zero-modes 



Here, we review the fermion zero- modes, which satisfy Eg. (12.31) with m„ = [IJ. On T 
the spinor ipn has two components. 
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= (ty) ■ (2.11) 



We use the gamma matrices on T as 



r- = ( , i ) . = ° S] . (2.12) 



1 Oy/' \i 

Then, the zero-mode equation is written as 

= 0, D^^.^o = 0, (2.13) 

where 

D = — (d + q-^{z + 0), (2.14) 



nR V ^ 2Imr 
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for the spinor with U{1) charge q. The charge q and magnetic flux m should satisfy that 
qm = integer. They also satisfy the following boundary conditions, 



(2.15) 



When qm > 0, only the zero-mode V'+.o has a solution, but ip-^ has no solution. Then, 
the chiral spectrum for the zero-modes is realized and the number of zero-modes is equal 
to qm. Their zero-mode wavefunctions are written explicitly as 



"(z + C) 



(2.16) 



where 



ef"^(^ + C,r)=exp 



'lm{z + C) , j , 

Trgmlmr 1 h 

Imr qm 



(2.17) 



+ i7iqmRe{z + Q [ ^"^[^ ^ + 2 ( ^ + 

Imr \ qm 



+ inqmHer { — — h 
qm 



Note that the effect of the Wilson line ( is the shift of the wavefunctions ip^''^™'{z) to 
il)h'i™'[^z + Q. The zero- mode wavefunction can be written by a product of the Gaussian 
function and the Jacobi ^^-function, i.e., 



"^(z + C) 



j/qm 



Z7l 



qm{z + ()lm{z + () 



Imr 



{qm{z + C),gmr), 



where 



,t) = ^exp [-Kiia + efr + 27ri{a + + b)] . 



(2.18) 



(2.19) 



('?/'+''™')* represents the anti-particle of '^+'''"5 ^ind is obtained from Eq. fl2.16p by replacing 
Ql''^"^{z + (, t) with Qj-''~'^"^{z + (, f). These zero-mode wavefunctions satisfy the following 
orthonormal condition. 



— I j,qm / I k,qm\ 



rp2 



(2.20) 



When qm < 0, there appear the zero-modes for ip-^, but not for ip+fi- The number of 
their zero-modes is equal to \qm\, and their wavefunctions are obtained similarly. In the 
following discussions, we assume qm > 0. 
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2.1.2 Fermion massive modes 

Here, we study the fermion massive modes witli m„ 7^ in Eq. f l2.3p . For m„ 7^ 0, the 
zero-modes, ip+n and ip-n-, mix each other in Eq. (12. 3p . They satisfy 

The 2D Laplace operator is defined as 



/\ = ]^{D\D}, (2.22) 
and it satisfies the following algebraic relations, 

|A.Z).] = ±^Dt = (2.23) 

Thus, massive modes are eigenfunctions of the Laplace operator A, and their mass spec- 
trum is derived in an analysis similar to the quantum harmonic oscillator. It is convenient 
to use the normalized creation and annihilation operators. 



^ A a^ = J-^D\ (2.24) 



ATiqm y Anqm 

which satisfy [a, a^] = 1. Then, the eigenvalues of the Laplace operator A are given as 

A„ = 27r^(2n + l), (2.25) 
and eigenvalues are also written as 

ml = An^n. (2.26) 
The corresponding wavefunctions ipn are written by 

^ ' J\njJ.qm 



Explicitly, the wavefunctions of massive modes are written as 



a^)"V'+7- (2-27) 



_ (2mgImr)V^ ^ 
^" ~ (2"ra!^)i/2 2^^^ + 



Im(2; + C) , j 



x/7„ V27rgmImr ; + £ 1 (2.28) 

\ \ Imr qm J J 

where Hn{x) is the Hermite function. Massive spectra of V'+.n and ip-^n are the same 
and the number of each of them is equal to qm. Note that ipif"^ satisfy the boundary 
conditions (I2.15p . Also, these wavefunctions satisfy the following orthonormal conditions, 

/ (i^(i#^'«™(^^''"^)* = (2.29) 

JT2 
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2.1.3 Scalar and vector modes 

Here, we study the scalar and vector modes on T^. The scalar fields are expanded as 
eigenfunctions of the Laplace operator, 

A<Pn{z) = ml<Pn{z). (2.30) 

That is, these eigenvalues are obtained as A„, i.e., 

= A„, = 27r^(2n + l), (2.31) 

for the scalar field. All of them including the lightest mode with n = are massive. 
Eigenfunctions are the same as those for the fermion, i.e. V'n'^'" Eq.f l2.28|) . 

Next, we study the vector field on T^. We are interested in the charged vector field 
with the U{1) charge g, where q ^ o|^ For example, they correspond to the vector 
bosons in the SU{2) gauge theory. We decompose the vector fields as Eq. f l2.ip . From 
Eq. (IB. lip in Appendix [Bl the mass-squared matrix is written by 

in the real basis of the 2D vector field, (0n,i,0n,2)- Instead of the real basis, we use the 
complex basis. 



1 X 1 

^(0n,l + «0n,2), 0n,f = "^v 



'^n,z = ^(0n,l + «0n,2), 0n,f = ^(0n,l " ^0n,2)- (2.33) 



The mass spectra of these internal wavefunctions are obtained through solving the fol- 
lowing equations, 

A - ^n,z = ml(t)n,z, (2.34) 

A + ) K-z = mlK-z- (2-35) 



That is, the mass spectrum of 0„^2 is obtained as 



r, , Anqm am , ^ , , 

ml = Xn- — |- = 27r^(2n - 1), (2.36) 



while the mass spectrum of is obtained as 



ml = X^ + ^ = 2vr^(2(n + 1) + 1). (2.37) 

The spectrum of (f)n,z includes the tachyonic mode for n = 0, while all modes of (j)n,z are 
massive. Their wavefunctions, (f)n,z and (f)n,z, are the same as those for the fermion, i.e. 
^mni Eq.([22HD- 



Obviously, there is no effect due to magnetic fluxes in the neutral vector fields with q = 0. 
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2.1.4 Massive modes only due to Wilson lines 

The massive modes also appear only due to non- vanishing Wilson lines a without magnetic 
fluxes. For completeness, we show their mass spectrum and wavef unctions. The internal 
wavefunctions for the spinor field as well as the scalar and vector fields satisfy the same 
boundary condition as Eq. fl2.15p with 



Xi 



Ima, 

Imr 



TT 



X2 



Imr 



-Imra. 



Then, the wavefunctions satisfying this boundary condition are obtained as 
1 



m ( 



exp 



/Ima \ Im^; 

iTx I V 2nji J Re^; + in- — - (— Rea + 2{nj — Rern/j)) 



\lmr Imr 
where ur, rij are integers. Their masses are given as 
o 47r^lmr 



m 



A 



+ 



1 



Imr 



/ Ima ^ ^ 
\Imr ^ 

2 

(— Rea + (n/ — Rernfi))' 



(2.38) 



,(2.39) 



(2.40) 



2.2 

Here we study the field theory on (T^)^. It is straightforward to extend the analyses on 
and (T^)^ to one on (T^)^. We use the complex basis, z"- = y^*"^ + r*2/^* with z = 1, 2, 3 
on the i-th T^, and the metric is written by 



ds^ = ^2{2ti R'fdz'dzK 



(2.41) 



We identify the complex coordinate as ~ + 1 and 2;* ~ 2;* + r*, and the area on the 
i-th is written by A' = An^ {R'^lmr'. 

We introduce the U(l) magnetic flux on the i-th as 



m ■ 
Imr' ' 



(2.42) 



where gm* is integer. This magnetic flux is derived from the following vector potential. 



A, 



IT 



-m^z\ 



A, 



TT 



2Imri"^ ' 2Imr 
We also introduce the Wilson line on the i-th T^, 



Tm^z\ 



(2.43) 



(2.44) 



Obviously, the mass spectrum and wavefunctions on each are given as those in 
section 12.11 The full eigenfunctions are the products of the eigenfunctions for the nj-th 
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modes on the i-th T^, and the full mass squared is the sum of masses squared for each 

. The number of massless fermions are obtained as n^gm*. The scalar field on 
is always massive. The vector field (p^r along the r-th (complex) direction on has the 
lowest mass squared with n* = (z = 1, 2, 3) as 



; \ ^ m m \ , ^ 



For example, when m^/A'^ + m^/A'^ — m^/A^ = 0, the massless mode appears in (p^ 
When / A? + / — m^/A^ is positive (negative), it becomes massive (tachyonic). 



3 Couplings including massive modes 

Here, we study couplings including zero-modes and higher modes in 4D low-energy effec- 
tive field theory. The 3-point couplings among zero- modes are computed in [1] |20] , and 
higher order couplings among zero-modes are studied in [T3j. First we briefly review on 
them in section 13.11 , and extend to 3-point and higher order couplings including higher 
modes in section 13.21 In section 13. 3[ we also consider the couplings including massive 
modes due to only the Wilson line effect, but not magnetic fluxes. 



3.1 Couplings among zero- modes 

Here we concentrate on the theory. We consider the coupling among three zero-modes, 
whose wave functions are given as (^ -)_ (i,t), ■?^*'<?2™'2(2; + (2,'^) and (?/'*''^^'"^(z -|- 

^3,r))*. They have U{1) charges, qi, q2 and q^, respectively, and the magnetic fluxes, 
mi, 1712 and appear in their zero- mode equations. We use the notation, A^^i = qimi, 
N2 = q2^2 and = qsin^. We assume that Ni,N2,N3 7^ 0. The gauge invariance 
requires that gi + ^2 = qs, Ni + N2 = N3 and Ni(i + N2C2 = N^^s. Their 3-point 
coupling in the 4D low-energy effective field theory is given by the following integral of 
wavefunctions, 

up to the 3-point coupling constant in higher dimensional field theory. Hereafter, we con- 
centrate on the part given as the overlap integral of wavefunctions, omitting the coupling 
constants in higher dimensions. For the Yukawa coupling, two of these modes correspond 
to the spinor fields, and the other corresponds to the 4D scalar field. The 4D scalar may 
be originated from the higher dimensional vector, e.g. on T^, if the 4D scalar is massless. 
At any rate, the wavefunctions are the same among the spinor, scalar and vector fields. 
Thus, we compute the 3-point and higher order couplings without specifying such Lorentz 
transformation behaviors. However, note that the Lorentz invariance leads to a certain 
selection rule. 
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In the computation of the above integral, the important property of zero-mode wave- 
functions is that they satisfy the following relation, 



N1+N2 



where 



X ^N2i-Nij+NiN2m, NiN2(Ni+N2) 



NiZi + N2Z2 Zi - Z2 

A — — — — , Y 



NI + N2 



(3.2) 
(3.3) 



N1 + N2 ' 
as shown in Appendix [C] (see also [H [20]). 

For example, when all of Wilson lines vanish, i.e. zi = Z2 = z, the above expansion 
becomes 



3,N2 I 



\A\Ni + N2) 



m=l 



i+j+Nim,Ni+N2, 



Z,T] 



X ^ 



■ N2i-Nij+NiN2m ' 
NiN2{Ni+N2) 





{0,tN,N2{N^ + N2)). 



(3.4) 



Then, by using the orthonormal condition fl2.20p . the 3-point coupling is obtained as 



_ f 2lmTNiN2 



+ N2) 



1/4 N1+N2 



,i+j+Nim 



m=l 



N2i-Nij+NiN2m ' 
NiN2{Ni+N2) 





iO,TN^N2iN^ + N2)). 



There is the selection rule for allowed couplings as 

k = i + j, (mod Ni] 



(3.5) 



(3.6) 



Similarly, we can calculate the 3-point coupling for non-vanishing Wilson lines. Its 
result leads to the 3-point couplings. 



y 



ijk 



+ N2 



+j+Nim^ 



ITT (iVi Ci ImCi +-;V2 C2lmC2 - A^3 Ca ImCs ) /Imr 



m=l 



N2i-Nij+NiN2m ' 
NiN2{Ni+N2) 





(3.7) 



(3.8) 



where the gauge invariance requires Ni + N2 + = A^4. For simplicity we consider the 
case that all of Wilson lines vanish, but it is straightforward to extend to the case with 



(iViiV2(Ci - C2), tN,N2{Ni + N2)). 
Next, we consider the 4-point couplings. 
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non-vanishing Wilson lines. The direct computation is possible by using the relation (jS 
However, the following calculation is much simpler [2]. We write the above integral 



= J dhd^z'ij''^'{z)ip^'''^{z)6\z - z')ilj'''^'{z'){t/j''''^{z'))y, (3.9) 

We replace the S function by 

6\z - z') = (3.10) 

s,n 

which is the summation over the complete set corresponding to eigenfunctions for the 
magnetic flux A^^i + This summation includes higher modes, i.e. n ^ 0. Then we can 
write 

s,n ^"^ 

x( [ rf2/^^.^l+^2(^')^'=.^3(^/)(^^,iV4(^/)))* j , (3 



Using the 3-point coupling among zero-modes, the above integral can be obtained as 

yijki ^J2y^^Y^~\ (3.12) 

s 

The higher modes n 7^ do not appear in this summation, because only zero-mode modes 
n = appear in the RHS of Eq.f l3.4p . 

Instead of Eq.f l3.9p . there is another way to split the integral, e.g.. 



y^Jke ^ j dhdh'^''^'{z)^^'^%z)6\z - z')^'''^^{z'){ilj^'^'{z'))y. (3.13) 
Then, by replacing the 6 function by 

5'{z - z') = 5:«'^^+^^(.))*<'^^+^^(.'), (3.14) 



t,n 

we obtain 



y'^''^ =J2y^''Y'^. (3.15) 

t 

We can show that both of Eqs.f l3.12|) and fl3.15p lead to the same result 
Similarly, we can calculate another type of the 4-point coupling. 



where the gauge invariance requires Ni + N2 = N3 + N4. Furthermore, the integrals for 
5-point and higher order couplings can be carried out in a similar analysis, and they are 
written by the proper summations over products of 3-point couplings. 
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3.2 Couplings including higher modes 

Here, we study couplings including higher modes. The relation f l3.2p among zero-mode 
wavefunctions plays an important role in the computation of the 3-point couplings for 
zero-modes. When we operate (9^, - |^^i)"H^^2 - i^^s)"' on the LHS of Eq. (Q, 
we obtain 



2Imr / \ 2Imr 

«f(^..-)-«f(^=.r). (3.17) 



On the other hand, when we operate {dz^ — ^^^i)'^^{(^z2 ~ ii^-^^a)"^ on the RHS of 
Eq. (13.21) . we obtain 

_ _L_iV,iV2(iVi + N^)y\ ^m-N,,+N.N2m,N,N2iN,+N2)^y^ 
, i+i+iVim,Ari+Af2/ V jN2i-Nij+NiN2m,NiN2{Ni+N2)(^ x /q i o\ 

^ ri+s \^}^) ' rni+n2~e-s yo.iQ) 

by using the derivatives with respect of X and Y. By identifying Eqs. ( 13.17^ and (I3.18p . 
the product of higher modes, ^/'^'^^(^i, r) and il>ii^'^{z2, r), is expanded acl 

A^l+Afa ni n2 ^{n2+i-s)/2 ^{ni-i+s)/2 



m=l 1=0 s=0 



{I + S)\{ni + n2 - I - s)\ ii+j+j^^rn,Ni+N2^v \ > N2i~Nij+NiN2m,NiN2{Ni+N2) n 

(3.19) 

■^A similar relation has been derived in twisted tori [16J. 
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When we take zi = z + (i and 2:2 = z + C2, it is found that 
V'^f (^ + Ci,r)-<^^(^ + C2,r) 



- \^ p r ( n»2-s ^i ^2 {£ + sy.jni + n2 - i - s)\ 

m=l 1=0 s=0 ' 

^Vi+s 12; + c,3, rj ■ ^ N2 [o.^O) 

Note that the last factor, ^^'"'^^^^(^^+^^)(^|^, r), is constant. 

Using the above relation, we can compute the 3-point coupling. 



Cn.na = y rf^rf#r(^ + Ci,r) -^^f (^ + C2,r) ■ «f^(z + C3,r))*, (3.21) 

in a way to similar to the 3-point coupling among the zero-modes. The result is obtained 

as 



N^+N2 ni n2 / ^^a+f-s Arm i , „ 

..ijfc _ Y^Y^ n r( ^\n2-s I ^^1 ^^2 [t: + s)l[ni +n2 - 1 - s)\ 

ynmana ~ / . / ,ni^e n2^s{ — i-) 



m=l e=o s=0 
> N2i-Nij+NiN2m, NiN2(Ni+N2) / Cl ~ C2 



(A^l + iV2)"i+"2+l ni\n2\ 



, iV2J — iViJ-|-iViiV2m, iViiV2l,iVH-JV2; / ^1 \ r r 



min(ni ,123) 
£=max{0,n:j —712 ) 



/^n2-n3+2^^„i+n3-2^ ^ _ 



JVafc-TVai, TVi jV2(Afi+jV2) ^ Cl - C2 x oox 

There is the selection rule among z, j and k, which is the same as one for the zero-modes 
(13. 6p . Thus, the flavor symmetry appearing only in zero- modes is still exact even when 
we take into account the effects due to higher modes. In addition, the following relation, 

n3 <ni+ n2, (3.23) 

should be satisfied for the mode numbers, ni,n2 and n^. For example, two zero-modes, 
rii = 77-2 = 0, can couple with only the zero mode = 0. On the other hand, the two 
zero-modes, 111=12^ = 0, can couple with higher modes, 77-2 7^ 0, and its coupling is 
determined by 



''^ j,N2k-N:ij, NiN2N3(_(l ^ t-\ 9A\ 

V7=h + N-2\N^ + N2) "^"^ Ui + iY,'"^^- ^^-^^^ 



Similarly, we can compute the 4-point coupling, 

yiiS2nsn, = I dzdz^ll^^iz + Cl, r) ■ <^^(^ + C2, r) ■ (z + Cs, r) ■ «f (z + (4, r))*. 

(3.25) 
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We rewrite it as 

d\d\'i^^^^{z + Ci, r) . i,i;^-{z + C2, r)5\z - + Cs, r) ■ (V^^f + (4, r))\ 

(3.26) 

and replace the 5 function by 

5'{z -z') = J2 rn'^'^'^'i^ + C, r)(C'^^+^^(^' + C, r)Y. (3.27) 

n,s 

Then, the 4-point couphng is given as the summation over products of 3-point couphngs, 

ijki ^ ijs ski /o 

i'nin2n3n4 / ^ yn\n2n ynn^rn' y^j.^Kjj 

n,s 

Similarly, we can compute other higher order couplings by products of the 3-point cou- 
plings. 

3.3 Couplings including massive modes only due to Wilson lines 

In the previous section, we have considered the couplings including higher modes under the 
magnetic flux. Here, we consider the couplings including massive modes only due to the 
Wilson line. The wavefunctions of such modes are obtained in Eq. (I2.39p . We compute the 
following 3-point couplings among two zero-modes ip^'^'^{z + Ci, t) and {ip'''^'^{z + (2, t))* 
and the massive mode ipn^^ni{z), where the two zero- modes have non- vanishing magnetic 
flux, while the massive mode has no magnetic flux, but Wilson line. Here, the gauge 
invariance requires that A^i = N2 and the Wilson line of the massive mode ipli^^niiz) 
satisfies A^iCi + 0^3 = -^iC2- Then, the 3-point coupling among these modes is given by 
the following integral, 

<).,n, = / dzd-zr'^^iz + Ci, r){^''''^{z + C2, r)r^^ZlM- (3-29) 
More explicitly, the integral is written by 



j,k 
y{W)niinj 



V2A^iImr 



2 



/ dzdz exp 



ITT 

{{Niz + ai)lm{Niz + ai] 



A^ilmr 



{NiZ + a2)lm{NiZ + a2)} + m(^-^ + ?j N^t -m(^-^ + n^ N-^f 

27rz|(iViZ + ai)(^^ + £^ - (A^i^ + ^2) + ^ 

f /Imaa \ Im^; / 
ni< Rez\ — h 2n/j 1 + ^ I — Reas + 2{ni — n^Rer) 



)) 



(3.30) 
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The integral over Kez imposes j = k — tir and i = n. In addition, the integral over Imz 
is Gaussian-like. By lengthy computation, it is found that 



-{(Imas + ri/jlmr)^ + (Reas — nj + n^jRer)^)} 



2iViImr 

(3.31) 



TT / _ 

+ i——-. Ima2Q;i + n/jlmrfai + ^2) + nnnilmr — n/Im(ai + 0:2 

Nilmr \ 

This behaves as a Gaussian function for the Wilson line a^. Thus, this coupling is sup- 
pressed depending on the Wilson line as well as and nj. The mode with the strongest 
coupling |z/({y)„^„J corresponds to the mode with npt = rij = 0, when 

1 Imas 11^1 , , 
- - < - — - < -, - - < Reas < -. 3.32 

2 - Imr - 2' 2 " ^ " 2 ^ ^ 

For other values of 0:3, another mode with non- vanishing riji and/or rij would have the 
strongest coupling. For example, for n/j = n/ = 0, we have 

|2 



\y(W)nR=nj=o\ - ^^^P 



2A^iImr 



(3.33) 



This coupling is suppressed depending on |a;3p. For example, when |a3p/2A'^iImr = 1, 
we obtain \y'i^^)nR=ni=o\ ~ ^"'^ ~ 0.04. The couplings to other modes with nR,ni 7^ are 
much more suppressed for the value of 03, which satisfy Eq. f l3.32p . 

Similarly, we can compute the 3-point coupling among two higher modes, ip-i^^^{z + 

Ci, t) and {ipn'j^^i^ + C2, t))* and the massive mode iplt^^miz), where the first two modes 
have non-vanishing magnetic flux, while the last mode has no magnetic flux, but Wilson 
line. We assume that ni < ^2. Such a coupling is obtained as 



y'nL(w)nRn, = J dzdz^^li^^ {z + (i, r)«^H^ + C2, r)WZlr (3-34) 

Again, the integral over Re^ imposes j = k — nn and i = n. For the integral over Kez, 
we use Eq. (IA.6p . Then, the result is written by 

(3.35) 



j,k y{w)nnni Y^^fc {riilf/^ f I / nn lma3 \Y 

yn,n,mnun, ^2^^7+7^^^ ^ - k)\ ^ V ^VTiViimT + ^^^^^ ^ ^ 



These couplings include the same suppression factor, yfw)nRni- 

Higher order couplings can be computed similarly. When we consider higher order 
couplings including more modes such as ipli^{T,z + we use the technique such as 
Eqs.( l3.26|) and fl3.27p . When we consider higher order coupling including more modes such 
as ip^^miz), we use the property that the product of two wavefunctions (-2) and 

iJ^mR,mi{z) is obtained as V'S^^.n^+mjl^), and that the Wilson line of V^i'^j.^^_„^+^^(2;) 
is just the sum of two Wilson lines, which (-2) and ip'^mR.m.iiz) have. Using these, we 

can compute higher order couplings. 
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4 Phenomenological comments 



We have calculated the couplings among zero-modes and higher modes. They have var- 
ious important implications from the phenomenological viewpoints. Here, we give some 
comments. 

The first example is about the proton decay. For instance, the proton decay would 
happen through the heavy X boson in the SU{5) GUT model. It couples with quarks and 
leptons by the gauge coupling before the gauge symmetry breaking. This coupling does 
not change in the 4D GUT theory even after the SU{5) group is broken. However, it can 
change in extra dimensional models, which have been discussed so far. Let us consider 
the SU{5) X f/(l) GUT model with extra space, or T^. We introduce non-vanishing 
magnetic flux m along the extra U{1) direction. Suppose that the 5 matter field has a 
charge q under the extra U{1) symmetry. Before SU{5) breaking, both of the quark and 
lepton in 5 are quasi-localized at the same place, and their coupling to the X boson is 
given by the gauge coupling. Then, we assume non-vanishing Wilson line a along the 
U{1)y direction in SU{5). It breaks the SU{5) gauge symmetry, the X boson becomes 
massive and its profile is written by 'ip^Lni=oi^) Eq. (12.391) . The quark and lepton in 
5 are still massless, but their profiles split each other, because of Wilson lines. In this 
case, the coupling among the quark, lepton and the X heavy boson is not equal to the 
gauge coupling, but it includes the suppression factor, |l/(^)„^=„^=ol' computed in the 
previous section. That is important to avoid the fast proton decay. For example, when 
|a;3p/2A''iImr = 1, we obtain |2/((^)„^=„^=ol ~ e~'^ ~ 0.04. Similarly, the couphng of the 
X boson with quarks and leptons in the 10 matter field can be suppressed. Then, the 
proton life time would drastically change by (9(10^ — 10^). 

Similarly, we can study the case that SU{5) is broken by the magnetic flux along the 
U{1)y direction, and the X boson becomes massive due to the magnetic flux. In this 
case, the coupling of the X boson with quark and lepton has the suppression factor given 
by Eq.( Km . 

Let us comment on another example. The Higgs mode gains its mass by non-vanishing 
Wilson line in certain models. That is, the Higgs mode corresponds to the open string 
stretching between two parallel D-branes ( on at least one of (T^)^) in the picture of in- 
tersecting D-brane models (see e.g. [21] ). The Yukawa couplings between this Higgs field 
and massless matter fields include the factor, \y'l^^-^nR=ni=o\- When the compactification 
scale is high such as the GUT scale and the Planck scale, this Wilson line a generating 
the Higgs mass is quite small, a <C 1, and the factor W(w-jn^=ni=o\ ^(1)- ^^^o 
important to see the moduli dependence of these couplings, that is, the dependence of the 
complex structure r and the Wilson line, a, which is the open string modulus. If F terms 
of complex structure and/or Wilson line moduli are non- vanishing, the corresponding A 
terms would appear and they are determined by the moduli dependence. At any rate, it 
is important to have found that the explicit moduli-dependence of these couplings, even 
though their values are of C(l). 

We would need heavy right-handed neutrino masses for the see-saw mechanism. These 
masses may be generated by non-perturbative effects (see e.g. [221 123]). Alternatively, 
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the right-handed neutrino masses are generated by Wilson hnes. If such a mass scale is 
comparable to the compactification scale, the couplings of the right-handed neutrino with 
the left-handed neutrino and the Higgs scalar would be suppressed. 

Finally, we comment on the Kahler metric. The Kahler metric of the matter fields is 
diagonal in the flavor basis. However, they couple with massive modes. Such couplings 
may induce off-diagonal entries in Kahler metric after integrating out massive modes. 
Such off-diagonal entries may lead to large FCNCs in the gravity-mediated supersymmetry 
breaking scenario [19]. However, when those couplings among massles modes and massive 
modes are suppressed, off-diagonal entries would be small. We have shown that the 
selection rule for allowed couplings including higher modes is the same as the one for only 
zero-modes. Thus, if there is a non-Abelian discrete flavor symmetry in massless modes 
|151 IT6] . that forbids the off-diagonal entries in the Kahler metric. Recall that such a 
symmetry is not violated by effects due to massive modes. 

5 Conclusion and discussion 

We have studied the mass spectrum and wavefunctions of zero-modes and higher modes 
in extra dimensional models with magnetic fluxes and Wilson lines. Furthermore, we have 
computed 3-point couplings and higher order couplings included higher modes in the 4D 
low-energy effective field theory. These couplings have non-trivial behaviors, because 
wavefunctions of massless and massive modes are quite non-trivial. Using our results, 
we can write down the 4D low-energy effective field theory with the full modes. Higher 
modes do not violate the coupling selection rules among only zero-modes. Thus, the flavor 
symmetry for zero-modes remains exact even when we take into account the effects due 
to higher modes. 

Our results would be important to phenomenological aspects, where couplings between 
massless and massive modes play a role, for example, the proton decay, the Higgs mass 
term, right-handed majorana neutrino mass term, FCNCs, etc. We will study in detail 
phenomenological applications of our results elsewhere. Threshold corrections and their 
moduli dependence after integrating out the massive modes would be important. 
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A Hermite function 



Here we show properties of the Hermite function, Hn{x), which is defined as 



Its derivative satisfies 



-^Hn{x) = 2nHn-i{x). 



The orthonormal conditions are written as 

1 

Hn{x)H^{x)e-''^dx = (5^,„2"V7rn!. 
We compute the following integral, 

/oo 
dxHnix + A)Hm{x + S)e-^^+^+^)', 
-oo 

for n < m. This integral can be calculated as 
J ox" 

/n 
fc=0 

n „ 

= J]2^C, „P.(-2S)'^-^ / dxe-^^^^^^^' H^J^Ux B) 

k=0 

r Jm—k 

ik=0 ^ 
n „ 

= J]2^Cfe „Pfe(-25)"-'=(-2A)"*-*^ / dxe-(^+^+^)' 
fe=o 

n 



fe=0 



where 



= ix) C^ = - Pt = — 

The following integral along the proper path, 

/oo 
-CX3 

leads to the same result as the above. 
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B Vector field 

Here, we study the (4 + 2) dimensional U{N) non-Abelian gauge theory (see also [HE]). 
Its Lagrangian is given as 



where 



MN 



OmAn - On Am - i[AM, Ajy] 



(B.l) 



(B.2) 



We compactify the two dimensions on with magnetic fluxes along U{1) directions. We 
decompose the U{1) parts and off-diagonal parts Wm, 



with 



Am = Bm + Wm = BljU^ + W^^Ca,, 



(B.3) 



(B.4) 



where Wm = (Wm)*. The quadratic terms of Wm in the Lagrangian are relevant to our 
study, and these appear 

1 



C 

where 



2r 



Tr [DmWnD^W^ - DmWnD^W^ - iGMN[W^'^, W^]) + 



(B.5) 



Gmn = QmBn — OnBm, (B.6) 
DmWn = OmWn - t[BN, W^]. (B.7) 

Here, the ellipsis denotes irrelevant terms. Furthermore, these terms are written by 



4^2 V V ■'■31 \ I 2g^ 



2 ( biWl^D^'W''' 



DiW-^D^W 



i,ab 



+ 



where 



We expand 



D,Wf = diWf - i{B^ - B\)W^ 



Wf{x,y) = Y,<^{^)K^{y)■ 



,(B.8) 



(B.9) 



(B.IO) 



Then, by imposing the gauge-fixing condition D^W^^ = 0, the equation of motion in the 
internal space is written by 



ab,i\ 



,2 lab 



(B.ii: 
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C Theta function identities and products of zero- 
mode wavefunctions 



Here we study the product of theta functions and product of zero-mode wavefunctions 
(see also [HEO]). The theta function satisfies the following identity, 







N2 




m=l 



r+s+Nim ' 
N1+N2 





[Zi + Z2,t{Ni+N2)) 



X 



N2r-Nis+NiN2m ' 
NiN2{Ni+N2) 





iz^N2 - TN^N2iN^ + N2)). (C.l) 



Using this identity, we can derive the following relations among products of zero-mode 
wavefunctions, 

N1+N2 



J,N2 I 



1 



ViVi + N2 



m=l 



+j+Nim,Ni+N2 



NiZi + N2Z2 
N1 + N2 



, T 



y. ^N2i-Nij+NiN2m,NiN2{Ni+N2) 



f Zi - Z2 

\N^ + N2 



, T 



(C.2) 



Its proof is as follows. The LHS is explicitly written as 



j,N2 I 



exp 



X ^ 



V ^2 

(iViZi,iVir)-^9 



(A^i^ilm^i + N2Z2lmz2) 



Ni 





J 

N2 




Imr 

{N2Z2,N2t) 



(C.3) 



and it can be rewritten by use of Eq. (]C.1|) as 



/2Imrv/]Vi7V2 

N1+N2 



1/2 



exp 



{NiZilmzi + A^2 -221^1-22) 



X 



N1+N2 



m=l 

N2i-Nij+NiN2m 
X NiN2{Ni+N2) 





Imr 

{N1Z1 + N2Z2,t{Ni + N2)) 

iNiN2{zi - Z2), rN^N2{N^ + N2)). (C.4) 



Since the exponent part is written as 



Imr 



{NiZilmzi + N2Z2lmz2 



Imr 

+ NMzi - Z2)lm 



(iViZi + N2Z2)lm 
zi - Z2 



( 



Nizi + N2Z2 



V N1 + N2 



Ni + N2 



(C.5) 



The RHS in Eq. (1C.4P is the summation over the products of wavefunctions, 

^i+j+Nim,Ni+N2 ^ N1Z1+N2Z2 ^■j ^N2i-Nij+NiN2m,NiN2{Ni+N2) ^ 



Z1-Z2 
Ni+N2^ 



T . 
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